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ABSTRACT 

The association of long gamma-ray bursts with supernovae naturally suggests that the cosmic GRB rate 
should trace the star formation history. Finding otherwise would provide important clues concerning these 
rare, curious phenomena. Using a new estimate of Swift GRB energetics to construct a sample of 36 luminous 
GRBs with redshifts in the range z = — 4, we find evidence of enhanced evolution in the GRB rate, with 
^4 times as many GRBs observed at z«4 than expected from star formation measurements. This direct 
and empirical demonstration of needed additional evolution is a new result. It is consistent with theoretical 
expectations from metallicity effects, but other causes remain possible, and we consider them systematically. 
Subject headings: gamma-rays: bursts — cosmology: theory 
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1. INTRODUCTION 

Long gamma-ray bursts are perhaps the grandest spec- 
tacles in astrophysi cs. Their connection with core- 
collapse supernovae (IStanek et al.ll2003t iHjorth et al.ll2003h 
indicates progenitors that are very massive, short-lived 
stars, leading one to expect the cosmic G RB history to 
follow the star formation history (SFH) ('Totani' Il997t 
Wiiers et a l.ll998l:lLa mb & Reichart 2000; Blain & Nataraja^ 
20001: iPorciani & Madauii2001 ) . Any deviation from this ex- 
pectation would provide new information about why a star 
should die as a GRB, complementing microphysical investi- 
gations (see lMeszarosI (1200 6) for a revie w). 

In just the past few years. Swift (iGehrels et al.l 120041) 
and a worldwide network of observers have detected 
gamma-ray bursts from higher redshifts than was pre- 
viously possible, spar king renewe d inter e st in the GRB 
redshift distributio n (|Berger et al.l l2005t iNataraian et al.l 
20051: iB romm & Loebl 120061 iDaigne et all l2006t 
Jakobsson et al. 2006; 'Le & Dermer"2007'; Yuksel & KistleI^ 
2007; Salvaterra & Chincarini 2007; Liang et al. 2007; 
Guetta & PiranI 120071: ICharv et al.1 12007,) . The Swift bursts 
with known redshifts now provide a sufficiently large sample 
upon which reasonable cuts can be made. Additionally, im- 
proved star formation measurements, compiled in the analysis 
of Hopkins & Beacom (2006), now provide a well-defined 
baseline for comparison. Using a simple, model-independent 
test, we find that the GRB rate history appears to evolve more 
strongly than the SFH until at least z w 4. We emphasize that 
our result establishes the existence of an evolutionary trend, 
but that, at present, we cannot identify the underlying reason 
(i.e., the distinction between correlation and causality). We 
discuss possible causes, including observational effects, the 
suggestion that the lower metallicity of the high-z universe 
could allow more bursts, and other possibilities that may 
result in additional GRB progenitors. 

2. THE EXPECTED GRB RATE 

An ever-increasing amount of data has brought about a 
clearer picture of the history of cosmic star formation. As 
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shown in Fig. [T] after a sharp rise up to z « 1, the SFH 
is nearly flat until z « 4 with relatively small uncertainties. 
These measurements are well-fit by a simple pi ecewise power 
law parametrization (iHopkins & Beacomll2006) . 
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scaled to psf(O) = 0.0197 Af© yr^^ Mpc"^ (a rate per co- 
moving volume). We parametrize the intrinsic (comoving) 
GRB rate relative to the SFH as riQ^Biz) — £{z) x psf(z). 
The fraction of bursts that can be seen at a given z, < F(z) < 
1, depends on the ability to detect the initial burst of gamma 
rays and to obtain a redshift from the optical afterglow. We 
cast the distribution of observable GRBs as 



dN ^ ^ £{z)psF{z)dV/dz 

dz (/beam) 1 + 2 



(1) 



This being an observed rate, cosmological time dilation 
requires the (1 + z)~^. The comoving volume ele- 
ment (in terms of the comoving distance, dc), dV/ dz = 
An (c/Hq) dl{z)/y^{l + z)3 + 17a, peaks at z 2.5. 
Dividing dV/ dz by the 1 + z term yields a volumetric fac- 
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Fig. 1. — Top panel: The history of star formation. Shown are the collection 
of data (circles) and fit of Hopkins & Beacom (2006). Bottom panel: The 
"volumetric factor" in Eq. (T). which distorts the observed GRB distribution 
when viewed in z. 
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Fig. 2. — The l uminosities (Ljso) o f 63 Swift gamma-ray bursts, determined 
from the data of IButler et"ai] 1200711 . vs. z (two 2:<0.1 GRBs are below 
10*** erg s~^). Points are plotted linearly in Q{z), as given in Eq. 0, to 
account for the volumetric factor. Open symbols may have underestimated 
Liso ■ The shaded region approximates the effective Swift detection threshold. 
Redshifts were measured as denoted, with medians of z = 0.8 for emission 
and z = 2.9 for absorption. 

tor that peaks at z ^ 1.4 (see Fig. [T]i"*. For a constant 
F{z), relatively fewer bursts should be observed at z ~ 4 
than z ~ 1. GRBs that are unobservable due to beaming 
are accounted for through (/beam) (see Bloom et al. (2003); 
iGuetta etal. (2005); Pan aitescul (120071) : iKocevski & Butled 
( 12007.) ; . Nava et al. (2007)). 

3. THE SWIFT OBSERVATIONS 

Swift has enabled observers to extend the reach of GRB ob- 
servations greatly compared to pie-Swift times, resulting in 
a rich data set. We consider bursts from the Swift archive^ 
up to 2007 May 15 with reliable redshifts and durations ex- 
ceeding Tgo > 2 sec. We estimate each GRB luminosity as 
Lisa = £'iso/[79o/(l + z)], where Ejso is the isotropic equiv- 
alent (beaming-uncorrected) 1 — 10* keV rest-frame energy 
release and Tgo is the interval observed to contain 90% of the 
prompt GRB emission. To form a uniform Li^o set, as dis- 
played in F' i g. [2l w e use the Eiso and Tgo values given by 
iButler et al] ( |2007|) . Note that Liso for several bursts (open 
symbols) may be underestimated due to Tgo values that over- 
estimate the GRB duration. 

To make the visual density of GRBs in Fig. |2] more mean- 
ingful, we define a linear x-coordinate that "corrects" for the 
effects of the volumetric factor as 



dz 



, dV/dz' 
l + z' 



(2) 



which we will quote in terms of (c/Hq)'^ « 79 Gpc'^. Remov- 
ing the influence of the volumetric factor allows for a bet- 
ter "by-eye" view, since dN/dQ = {dN / dz) / {dQ / dz) = 
F{Q) ?i-GRB(Q)/ (/bcam)^ SO that a flat GRB rate would ap- 
pear as a constant density of GRBs with similar Liso pei" linear 
interval in Q. 

The Swift trigger is quite complex, working in the 
15 — 150 keV band and using time-dependent background 
subtraction and variable time windows in order to maximize 
burst detection (Gehrels et al. 2004; Band 2 006). While t his 
sensitivity is "very difficult if not impossible" (lBandll2006l) to 

^ Using Ha = 0.7, Qm = 0. 3, and Hq = 70 km/s/Mpc; changing these 
requires correcting the SFH (see lHopkinsI 120041) ) 

' http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table 



Fig. 3. — The differential GRB distribution versus Q. Outlined bins contain 
the set of all 63 GRBs (median z = 2.3), while shaded bins contain the 44 
bursts with Liso > 10^^ erg s^^ (with a median of z = 2.9). The rise seen 
suggests that GRB rate evolves more strongly than star formation. 



parametrize exactly, an effective luminosity threshold appears 
to present in the data (roughly estimated as oc in Fig.|2]i. 

A representative sample of bursts in a given redshift range 
can be selected, while avoiding detailed assumptions concern- 
ing this threshold, by simply choosing a lower cutoff in Liso- 
Considering only bursts that could have been seen from any- 
where within the redshift range examined allows for the Swift 
contribution to the F{z) term to be treated as constant in 
z. This technique effectively integrates the GRB luminos- 
ity function, dN/dLiso, above the chosen Ljso cut without 
assuming its functional form. This reduces the problem to 
number counts. 

4. TESTING FOR EVOLUTION 

Because psf{z) is now reasonably well-measured from 
z « — 4 (and nearly flat for 1 < z < 4), we consider GRBs 
in this range for comparison. Using only bursts with Liso > 
10^^ erg s^^ creates a set of GRBs with a minimal expected 
loss of GRBs up to z < 4. Fig.[3]displays the differential dis- 
tribution of these GRBs versus Q. The bursts above this cut 
(shaded bins) can be compared to the set of all bursts (out- 
lined). As can clearly be seen, removing the population of 
lower luminosity bursts that are only observable at low z re- 
veals a distinct rise in the observed number of "bright" bursts. 
The drop at z > 4 is likely due to the Swift threshold and an 
overall drop in star formation. 

To compare the GRB data with the SFH expectation, we 
make use of Eq. ([T]), parameterizing the "effective evolution", 
for simplicity, as F{z) x £{z) oc (1 + z)". Since the lumi- 
nosity cut removes the influence of the Swift threshold, any 
z-dependence of F would be due to other potential observa- 
tional effects. We compare the predicted and observed cumu- 
lative GRB distributions in Fig.|4] This compares the relative 
trends, independent of the overall normalization and the value 
of (/beam)- A Kolmogorov-Smirnov test reveals that the SFH 
fit alone is incompatible at around the 95% level. 

Positive evolution results in clear improvement, strength- 
ening indications based upon differen t analyses and smalle r 
data sets (e.g., Daigne et al. (200^; iLe & Derme"^ (l2007h : 
lYiiksel & Kistleil (i2007i) ). We find that the K-S statistic is 
minimized for a= 1.5. To simulate dispersion in the data, 
we replace the Tgn values used to calc ulate Ljso with those 
reported by Swift dSakamoto et alj2007h and find similar con- 
clusions. In neither case is any significant correlation of lumi- 
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Fig. 4. — The cumulative distribution of the 36 Swift GRBs with Liso > 
10^^ erg in 2 = — 4 (steps), compared to the expectations from an 
effective evolution of (1 + z)^-^ with respect to psp{z) (dark solid); and the 
SFH alone (dashed), which is inconsistent with the data at ~ 95%. 

nosity with z found i n this range, disfavo ring just an evolving 
luminosity function dEfron & Petrosianl !T992). Although we 
have attempted to minimize the loss of GRBs due to the Swift 
threshold, we are more likely to be underestimating burst 
counts at the upper end of the z-range examined, which could 
imply even stronger evolution. 

5. IS THE EVOLUTION INTRINSIC? 

This evolution could result from several causes. The K-S 
test disfavors an interpretation as a statistical anomaly. While 
Psf{z) could just be mismeasured, the relatively-small uncer- 
tainties over the range in question suggests that this is not the 
likely origin. We will discuss in more detail whether it is due 
to some selection effect (i.e., a changing F{z)), or a real phys- 
ical effect related to changes in the number of progenitors (an 
evolving £{z)). 

It has been suggested that a hig her percentage of GRBs 
may g o undetected at low redshifts (lBloomll2003llFiore et all 
|2007|) We perform several simple diagnostics to determine 
whether the sample of Swift GRBs without redshifts contains 
an overwhelming number of bright, low-z bursts. To simplify 
analysis, we consider a sample of bur s ts tha t meet the "de- 
tectability" criteria of iJakobsson et al.l (12006 ) (in particular, 
low Galactic extinction and quick Swift localization), a set of 
50 GRBs with a confirmed z and 47 without. With this set, 
using the procedure of Sect. |4] the SFH expectation alone is 
still incompatible at around 95%, with a =1.5 minimizing the 
K-S statistic. 

On average, more distant bursts are expected to have lower 
observed gamma-ray fluxes, JF, estimated by dividing each 
15 - 350 keV fluence and Tgo from Butler et al. (2007). In- 
deed, GRBs at z >2 have {Tz>2) ^ 1-4 x IQ-'^ erg cm^^ 
s~^, compared to {J-'z<c2) ^ 2.6 x 10^^ erg cm~^ s~^ at 
z < 2. For GRBs without a redshift, the average flux is just 
{J^no z) ~ 0.7 X 10"'' erg cm"^ s~^ A two-sample K-S 
test between the z < 2 and z-less sets reveals that they are in- 
compatible at ^-^70%. Limiting the z < 2 set only to GRBs 
with Liso > 10^^ erg s^^, this increases to ~99%. While not 
conclusive, these results could be interpreted as most z-less 
bursts being either at high z, in which case our evolution may 
be underestimated, or at low z with lower intrinsic luminosi- 
ties, which may not survive our Liso cut. Additionally, we ex- 
amine the fractions of GRBs detected by Swift's UV-Optical 
Telescope in the Swift archive. Of bursts with a known z, 34 



were seen by UVOT with (z) = 2.2, while the 16 not seen 
have (z) ~ 3.0. Bursts lacking redshifts seem more consis- 
tent with the high-z set, with only 8/47 seen by UVOT. 

We also consider whether the fraction of observable bursts, 
F{z) = fz f Swift, is somehow increasing with z. While 
f Swift is difficult to quantify, our selection criteria disfavor 
incompleteness of our sample at low z. We focus upon the 
probability of determining a redshift for a given GRB, which 
we subdivide as = /e/a /ag /human- Perhaps the most ob- 
vious influence on this term is the fact that at low z, most red- 
shifts are determined by observing emission lines (from the 
host galaxy); while at higher redshifts, nearly every redshift is 
found through absorption lines in the afterglow spectrum (see 
Fig.|2]and the Swift archive). While this emission/absorption 
bias, /e/A' is not easy to quantify, it should not cause such 
an evolutionary effect, since most redshifts in our sample are 
found through absorption. 

The observability of an optical afterglow, /ag. might be ex- 
pected to be steeply f alling with re dshift; however, for a spec- 
trum oc t^" v^f^ (Sari et al. 1998), cosmological redshifting 
may allow the earliest (brightest) portion of the afterglow to 
be more visible (Ciardi & Loeb 2000). If the flatness of the 
obscuration-corrected SFH at moderate z arose from a steeply 
increasing uncorrected rate and a deceasing dust correction, 
the detectable afterglow fraction might increase with redshift; 
this does not appear to be the case (Schiminovich et al. 2005). 
A paucity of IR detections of hosts of "dark" GRBs a lso ar- 
gues against a significant obscured fraction (.Le Floc'h et al.l 
l2006l) . The decision of which telescopes are made available 
to observe GRBs, along with other such non-intrinsic proper- 
ties, can be folded into /human- None of these terms appear 
to be able to increase overall observability with z, disfavoring 
an origin of the trend in F{z). 

6. POTENTIAL SOURCES OF EVOLUTION 

We now investigate whether an evolving E{z) can explain 
the observed evolution. While it is now generally accepted 
that long GRBs arise from massive stars, the special con- 
ditions that are required for s uch an event are still in ques - 
tion. In the collapsar model (MacFadv en & Wooslevlll999l) . 
the collapse of a rapidly-rotating, massive stellar core to a 
black hole powers a jet that is seen as a GRB. Since ev- 
ery obse rved supernova coinciden t with a GRB has been of 
Type Ic (IWoosley & Bloomll2006h . the progenitor star should 
also have lost its outer envelope (without losing precious 
angular momentum). Rather strong observational evidence 
now indicates th at GRB host galaxies tend to be faint an d 
metal-poor (e.g. iFruchter et al) (12006"); jStanek et al.l (l2006h : 
iLe Floc'h et all (12003); Fynbo et al. (2001), increasing inter- 
est in models that use single, low-metallicity stars as a path- 
way to a collapsar ( Yoon & Langer 2005; Woosley & Hegei 
2006). Decreasing cosmic metallicity may cause the GRB/SN 
ratio to rise with z; e.g., the prediction of I Yoon et al.l (l2006h 
can be estimated as (1 + z)^ '^ dYiiksel & Kistled2007h . While 
a preference for low-metallicity environments may be the sim- 
plest explanation, absorption line metallicity studies remain 
inconclusive (iSavaglioi 120061 ; IProchaska et al.l 120071) . One 
may therefore wonder whether it is possible to concoct an evo- 
lutionary scenario without direct progenitor metallicity de- 
pendence. 

If GRBs are in stead p roduced in binary s ystem s 
(e.g.. iFrvereta l. (1999); ' Podsiadlowski et al.l (2004); 
iDale & DaviesT (2006)). some other mechanism might 
lead to the appearance of evolution. Since most GRBs 
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appear to occur in star clusters dPruchter et alj l2006l) . 
where the fraction of massive stars in binaries may be 
high (Kobulnickv & Frveil l2007h . such channels could be 
important. For example, the merger of two >15Mq stars 
with Mi/M2>0.95 (i.e., "tw ins", which may be common 
(iPinsonneault & Stan^l2006h ) in close orbits (r < few AU) 
can lead t o a more-massive, rap idly-rotating core lacking an 
envelope (iFrver & Hegerl l2"005b . In a dense cluster, close 
binaries tend to end up closer due to g ravitatio n al sca tterings 
with interloping stars (Heggie's Law; Heggi^ (Il975l) ). with 
a scattering timescale of < 10 Myr f or a stellar density of 
p 10'' Mq pc"3 (iHut et al.1 Il992h . An increased rate 
of "interloper-catalyzed" binary mergers (ICBMs) could 
result from a larger fraction of star formation occurring in 
such environments at higher z, and seen as an enhancement 
in the GRB rate. Such a speculative scenario presents 
a dynamical source of evolution (instead of altering the 
microphysics) and, as these rates are cx ^p^), possibly allow 
for examination of the "density contrast" of star formation. 
It is interesting that GRBs were discovered in searches for 
gam ma rays from explosion s related to ICBMs of a different 
sort (iKlebesadel et al.lll973l) . 

An evolving IMF, becoming increasingly top-heavy at 
larger z, would increase the relative number of massive stars 
produced. Since star formation measurements are primarily 
based on radiation from such stars, this alone may not lead to 
apparent evolution, unless the very massive end (> 25Afo) 
chan ged significantly. Any evidence of ev olution in the IMF 
(e.g., IWilkins et al.l (IMNRAS. submitted!) ') provides motiva- 



tion for considering such a change. However, this need only 
occur in those galaxies that host gamma-ray bursts. A sim- 
ilar effect could result from a larger high-z population of 
small galaxies resembling low-z GRB host galax ies, seen as 
evolution in the galaxy luminosity function (e.g.. lRyan et al.l 
(200 7)). An additional possibility is that, if the massive stellar 
binary fraction decreases with decUning meta lUcity, as might 
be evidenced in the LMC (iMazeh et al.ll2006h . the number of 
potential single-star progenitors may increase. 

In conclusion, the set of Swift gamma-ray bursts now al- 
lows for model-independent tests of the connection between 
the GRB and star formation rates. We present quantitative 
evidence that the GRB rate does not simply track star for- 
mation over a broad range in redshift; some mechanism, of 
a presently-unknown nature, is leading to an enhancement in 
the observed rate of high-z gamma-ray bursts. The effects 
of stellar metallicity appear to be a compelling explanation, 
but cannot be proven yet. Additional observations will be the 
only way to discern the root cause of this effect, and allow for 
proper understanding of the underlying astrophysics. 
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